Direct evidence for inversion formula in multifractal financial volatility measure 
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The inversion formula for conservative multifractal measures was unveiled mathematically a 
decade ago, which is however not well tested in real complex systems. In this Letter, we pro- 
pose to verify the inversion formula using high-frequency turbulent financial data. We construct 
conservative volatility measure based on minutely S&P 500 index from 1982 to 1999 and its inverse 
measure of exit time. Both the direct and inverse measures exhibit nice multifractal nature, whose 
scaling ranges are not irrelevant. Empirical investigation shows that the inversion formula holds in 
financial markets. 
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In recent years, the concept of inverse statistics has at- 
tracted much attention in turbulence [l], Q as well as in 
financial markets 0] based on time series analysis. The 
direct structure function concerns with the statistical mo- 
ments of a physical quantity g measured over a distance 
s such that S q (s) — (g\\(s) q ). The multifractal nature 
of direct structure functions has been well documented 
in turbulence 0, |j| Hj], as well as in finance 0, H, Q, 
which is characterized by S q (s) ~ s^ 9 -* with a nonlinear 
scaling function C(<z)- I n contract, the inverse structure 
function is related to the exit distance, where the physi- 
cal quantity fluctuation exceeds a prescribed value, such 
that T p (g) — (s p (g)). One can intuitively expect that 
there is a power law scaling stating that T p {g) ~ g^ p \ 
where 4>(p) is also a nonlinear function. Furthermore, if 
s ~ g^' h , Schmitt has shown that there is an inversion 
formula between the two types of scaling exponents such 
that = — p and 4>(p) = —q 10]. A similar intuitive 
derivation for the inversion formula is given for Laplacian 
random walks fill ]. 

The power-law scaling in inverse structure function was 
observed in the signals of two dimensional turbulence 
[3 EH ,in the synthetic velocity data of the GOY shell 
model jl|, 03] , and in the temperature and longitudinal 
and transverse velocity data in grid-generated turbulence 
[To] ]. However, this scaling behavior was not observed 
in other three dimensional turbulent flows from different 
experiments @, [l6|, [TtJ • The inversion formula for direct 
and inverse structure functions is verified for synthetic 
turbulence data of shell models [l4] but not for wind- 
tunnel turbulence data, which cover a range of Reynolds 
numbers Re = 400 - 1000 [TJ]. 

It is argued that [181 ] , the absence of inversion formula 
between the scaling exponents of direct and inverse struc- 
ture functions is due to the facts that the velocity fluctu- 
ation is not a conservative quantity while a strict proof of 



the inversion formula was given for conservational mul- 
tifractal measures [l9|, [2(| • It is noteworthy pointing out 
that, the inversion formula given by Roux and Jensen 
is obtained based on a special case of conservative mea- 
sures, although they verified the inversion formula in the 
direct and inverse structure functions of shell models [3] . 
It is thus natural that Xu et al. proposed to test the in- 
version formula in the energy dissipation rate (a kind of 
conservative measure) rather than in the structure func- 
tions and they did found sound evidence in favor of the 
proof |U|. 

The inversion formula was theoretically established for 
both discontinuous and continuous multifractal measures 
by Riedi and Mandelbrot [l9, 20]. Let /i be a probability 
measure on [0, 1] whose integral function M(a) — /i([0, a]) 
is right-continuous and nondecreasing. Since the measure 
is self-similar, we have /i = X)T=o TO *A t ( u 'i~ 1 ('))i where 
wi 's are the similarity maps with scale contraction ratios 
Ti G (0, 1) and 2t=i m i ~ 1 with m i > 0. The multifrac- 
tal spectrum f(a) of measure \i can be obtained via the 
Legendre transform of r, which is defined by 



2^ n*?r, r = 1 



(1) 



The inverse measure of fj,* can be defined as follows, 
' inf{a : M(a) > b} if b < 1 



H* = M*(b) = 



1 



if b = 1 



(2) 



where M*(b) is the inverse function of M(a). Since /i 
is self-similar, its inverse measure fi* is also self-similar 
with ratios r* = m, and probabilities m* = fj, whose 
multifractal spectrum /*(«*) is the Legendre transform 
of 9, which is defined implicitly by 



i 



The inversion formula follows immediately that 



'Electronic address: wxzhou@ecust.edu.cn 



r(q) 
6(p) 



-P 

-q 



(3) 



(4) 



2 



Equivalently, we have 

T{q) = -6-\- q ) 

or 

0( P ) = -T-\- P ) . 



(5) 



(6) 



These two equivalent relations are testable. Follow- 
ing this line, the inversion formula was verified in fhSj 
with high- Reynolds turbulence data collected at the SI 
ONERA wind tunnel [||, which is however the only ev- 
idence. Due to the well documented analogues between 
turbulent flows and financial markets [2l[, in this letter, 
we propose to test the inversion formula in financial mar- 
kets using high-frequency historical data of the S&P 500 
index. Our data consist of 18-year minutely prices span- 
ning from 1 January, 1982 to 31 December, 1999 with a 
total of 1.7 million data points. The minutely return r(t) 
is calculated as follows, 



r(t)=ln[I(t)/I(t-l)] 



(7) 



where {I(t) : t = 1, • • • , T} is the time series of minutely 
S&P 500 index. 

We first construct the direct volatility measure and in- 
vestigate its multifractal nature. The absolute return is 
utilized as a proxy for volatility such that v(t) — \r(t)\. 
According to the partition function method for multi- 
fractal analysis, the series is firstly covered by N boxes 
with identical size s = T/N. The sizes of the boxes are 
chosen such that the number of boxes of each size is an 
integer to cover the whole time series. On each box, we 
construct the direct measure as 



Hn(s) 



£ «(*). 

t=(n-l)s+l 



(8) 



where V = J2t=i V W an( ^ n = 1, ■ ■ ■ , N . By construc- 
tion, this volatility measure /x is conservative. The pres- 
ence of multifractality in /i has been confirmed based on 
the multiplier method utilizing the same data set of the 
minutely S&P 500 index [22j . Alternatively, the volatility 
measure fi of Chinese stocks and indexes exhibits multi- 
fractal behavior based on the partition function approach 

M- 

For order q, the direct partition function Xq( s ) can be 
estimated using 



N 



Xq 



(9) 



71=1 



When /i -C 1 and q 3> 1, the estimation of the parti- 
tion function \ will be very difficult since the value is so 
small that it is "out of the memory" . To overcome this 
problem, we can calculate the logarithm of the partition 
function lnx 9 (s) rather than the partition function itself. 
A simple manipulation results in the following formula 



InXgO) = ln£ 



Mn(s) 



gin fj,„ 



(10) 



where /z max = max{/i„ : n = 1, • • ■ , N}. This trick ap- 
plies for the calculation of inverse partition functions as 
well. 

Figure [1] plots [Xq{ s )] 1 ^ q ^ 1 as a function of box size s 
for different orders. Sound power laws are observed for 
each partition function such that 



Xq(s) 



Q r(q) 



(11) 



in which the scaling range spans about three orders of 
magnitude. The scaling exponent r(q) can be estimated 
through a power-law fit to the data in the scaling range. 
We will see that r(q) is a nonlinear function, confirm- 
ing the presence of multifractality in the direct volatility 
measure. 
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FIG. 1: (color online.) Dependence of Xq( s 
tion of box size s for q — —3, q = —1, q 
q = 6. The curves have been translated vertically by a factor 
of 0.001, 0.01, 0.1, 10, and 100 in turn for better visibility. 
The solid lines are power-law fits in the scaling range. 



We now investigate the scaling behavior of inverse par- 
tition function of exit times. For each threshold Ad, a 
sequence of exit times Sj(Av) can be determined succes- 
sively from j = 1 to j = J by 



£ s k = inf ji : J v{t)dt > jAv j 



(12) 



where v(t) = v(t) for t € [t, t + 1). The inverse measure 
is defined as the normalized exit time 



ti*(Av)= Sj /T , 



(13) 



and the inverse partition function can be determined as 
follows 



x;(Av) = J2[^( A v)] p 



where 



1 

Av 



v{t)dt 



(14) 



(15) 



3 



Figure[2]shows the dependence of xZiAv) 1 ' on the 
threshold values Aw for different values of p. Power-law 
scaling can be observed 



X *Js) ~ Av e ^ . 



(16) 



where the scaling range covers about three orders of mag- 
nitude. The straight lines are the best fits to the data in 
the scaling range, whose slopes correspond to the expo- 
nents 9(p)/(p — 1). We note that the two scaling ranges 
(si,S2) and (Avi, Av 2 ) for direct and inverse partition 
functions are related by 



Aw 



s x v„ 



(17) 



where w moan = 2.69 x 10 4 . Specifically, we find that 
Avi w si x w moan and Av 2 ~ s 2 x w moan . This puts 
forward sound evidence upon the determination of the 
scaling laws. 
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FIG. 2: (color online.) Dependence of Xp(Av) 1/,p_1 on the 
thresholds An. The same translated processes are also per- 
formed on the data points of p — —3, p = — 1, p — 0, p = 4, 
and p = 6. The solid lines are the best fits to the data. 

A subtle issue concerning negative moments arises, 
which is related to the probability density functions 
(PDFs) of volatility and exit time respectively. The 
volatility of S&P 500 index is log-normally distributed in 
the center followed by a power-law tail for larg e volatil- 
ities while the right tail seems truncated [24l [2o| . In 
addition, due to the construction of the volatility mea- 
sure, /i(s) > for large s in the scaling range shown in 
Fig. [2 Therefore, negative moments can be estimated 
numerically. Taking into account the statistical signifi- 
cance of the estimation of partition functions fl7l.l26l], we 
focus on q <E [—4, 8]. 

The PDFs of exit time defined in this letter have not 
been investigated before. Let us denote by f(s) the PDF 
of exit times for a fixed threshold Av. For comparison, 
we normalize the exit time by their standard deviation 
a(Av) for each given Av. Then the PDF of the normal- 
ized exit times x — sja can be determined by 



Figure[3]shows the empirical PDFs af(xa) of the normal- 
ized exit times x = sja for different thresholds Av. The 
PDFs at different thresholds cannot be superposed by the 
simple normalization procedure. As sketched in Fig.[3Ja), 
the PDFs are strongly asymmetric and not log-normal, 
which differs remarkably from the situation of energy dis- 
sipation in three-dimensional fully developed turbulence 
showing roughly log- normal distribution [181 ] . More inter- 
estingly, the probability density functions show plateaus 
on the left tails, which ensures the existence of any neg- 
ative moments. Figure [31(b) shows that the right tail 
relaxes exponentially for small thresholds or faster for 
large thresholds. This relaxation behavior is different 
from those exit times extracted from financial return se- 
ries exhibiting a power-law tail % S3, SI SI, SS S3 , 
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p(x) = af(xa) . 



(18) 
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FIG. 3: (color online.) Empirical probability density func- 
tion p(x) of the normalized exit time x = s/a for different 
thresholds Av = tWan/2 (o), u mcan (□), and 2« mcan (o). 

The power-law exponents r(q) for direct partition func- 
tions are plotted as open circles in Fig. [U while the expo- 
nents 9{q) are illustrated as triangles. Both r(q) and 9(p) 
are nonlinear, indicating that the time series of volatility 
and exit time possess multifractal properties. The func- 
tion — 6~ l (— q) is determined numerically from the 9{p) 
curve, which is plotted in Fig. 2] as a dashed line. We can 
find that the values of — 9~ 1 (— q) are in excellent agree- 
ment with the values of t(j>), which provides strong evi- 
dence supporting the inversion formula in Eq. Simi- 
larly, the — t -1 (— p) curve numerically obtained from the 
t(jp) function is depicted as a solid line, which coincides 



4 



remarkably with the 9{q) curve. In other words, the in- 
version formula Eq. ^ also holds as expected. We note 
that the differences between the comparing curves are 
well within the error bars. 
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FIG. 4: Testing the inversion formula in financial volatility. 

In summary, we have attempted to test the inver- 
sion formula for conservative multifractal measures us- 
ing high-frequency volatility data of the S&P 500 in- 
dex. We have performed multifractal analysis on both 



volatility and exit time series based on the partition func- 
tion method. Our investigation confirms that both direct 
and inverse partition functions exhibit nice multifractal 
properties. The two scaling ranges are consistent with 
each other. Furthermore, we found that the function 
— r~ 1 (— p) extracted numerically from r(g) overlaps with 
the 9{p) curve and the function —6 {—q) determined 
from the 9(p) curve collapses on the r(q) curve, which 
verifies the inversion formula. We also investigated for 
the first time the empirical distributions of exit time of 
financial volatility at different thresholds. The PDFs of 
exit time are nontrivial, which are neither log-normal nor 
power laws observed in other systems. 
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